Non linear transport properties of an insulating YBCO nano— bridge 
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We have investigated the transport properties of an insulating sub-micrometric YBa2Cu307_4 
bridge, patterned on a thin film. As expected for a variable-range-hopping insulator, transport is 
found non linear. The reduced dimension allows for the observation of an individual fluctuator 
generating random telegraph noise, which dynamics could be explored as a function of the tempera- 
ture and transport current. Some recordings clearly exhibit three levels fluctuating resistance, with 
comparable level separation and correlated dynamics, which likely result from the existence of two 
states or correlated clustered charge traps. 

PACS numbers: 71.30.-|-h,72.20.Ht,72.70.-|-m,74.72.Bk 



It is well known that superconductivity in cuprates 
emerge from doping a Mott insulator to obtain the under- 
doped superconductor As such, as pointed out early, 
high-Tc superconductors share similarities with doped 
two-dimensional weakly compensated semiconductors 0, 
0]. In the insulator, carriers are localized within a lo- 
calization length (effective Bohr radius) which diverges 
- together with dielectric constant - at the transition to 
the delocalized state, at some critical value of the tuning 
parameter, x = Xc, where x may be a carrier density, dis- 
order strength, impurity concentration, magnetic field... 
In the perspective of a quantum phase transition between 
adjacent insulating and superconducting states, the vari- 
able range hopping (VRH) law for the insulator conduc- 
tivity may be viewed as one example of a scaling law 
associated to a quantum critical point, in which the di- 
vergence of the localization length at the critical density 
for doping accounts for the critical temperature scaling of 
the form p{T,n) = F[5/T^/^''), where 5 = {x-Xc)/xc, z 
is the dynamical exponent and v is the correlation length 
exponent. In the regime of hopping in the Coulomb 
gap 01 1 ones has z = 1, as expected for a strongly inter- 
acting Coulomb system^. There are several experimen- 
tal indications for the classical metal-insulator transition 
that the VRH regime might be characteristic of the in- 
trinsic insulating phase, representing the disordered state 
to be found on one side of such a quantum transition, 
even in the absence of strong Coulomb interaction[^ Q. 
In the case of strongly disordered systems, the insulating 
phase could be found to coexist with the metallic or su- 
perconducting one [6, 8] and may influence the transport 
properties far into the metallic phase Besides classi- 
cal scaling dependence for the conductivity, the insulat- 
ing phase is also characterized by the existence of slow 
dynamics for charges localized on traps. In the modeliza- 
tion by a percolation resistance network, such a behav- 
ior arises from the existence of a network of thermally 
assisted hopping paths for conduction, and neighboring 
sites loosely coupled to this network pill Due to the 
exponential dependence of the hopping rate on site dis- 




FIG. 1: a) Schematics for the FIB process from a lithogra- 
phied fllm, b) actual realization on a YBa2Cu307-5 film, c) 
detail of the bridge. 



tance, there is a wide distribution of hopping rates from 
such sites, which accounts for a 1// resistance noise[l2!]. 
Reducing the dimension of the conductor may allow to 
resolve some individual contribution of the fluctuators to 
the resistance. In particular, sites with occupancy vary- 
ing slowly on the time scale of the experiment may be 
evidenced by random telegraph noise. 

We have investigated the transport properties of 
such a conductor, from a sub-micron bridge pat- 
terned on YBa2Cu307_5 thin film. There are indica- 
tions that cuprates high-Tc superconductors may pro- 
vide one example of a quantum critical point, at the 
superconductor-insulating transition [T3|. For weakly 
disordered compounds, such as the present one, there 
should be no mixing of the two phases, nor existence 
of some intermediate glassy phase at the transition. It 
is however, notoriously difficult to prepare such bridges, 
with a controlled doping state. The reason for this is 
likely the ability of the doping oxygen to migrate from 
narrow bridges, yielding patterned samples that are less 
doped than the starting material. Complex prepara- 
tion of the bridges, combining appropriate deposition 
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FIG. 2: Resistance of the constriction {1 = 1 nA), in repre- 
sentations where a VRH would appear linear (2D : T'^^'^; 3D 
: T"^/*). Both VRH laws fail below T « 45 K. 



of protective layers and reactive ion etching, such as in 
ref. p!3 |. could provide bridges with unchanged doping 
state. However, more direct techniques, such as focused 
ion beam (FIB) etching, likely induce a strong reduction 
of doping, so that the observation of superconductor- 
insulator transition with the reduction of the bridge di- 
mension may actually be due to the concomittent reduc- 
tion in the doping state|15]. The migration of the oxygen 
out of the material could be due in this case to both the 
increased local temperature during etching, and to elec- 
trical field effects created by the charged ions deposited 
by FIB. The film was grown on a SrTiOs substrate, us- 
ing laser ablation [l^. and pre-patterned using optical 
lithography. The as-grown sample displayed sharp su- 
perconducting transition (0.5 K wide) at T = 92 K; 
the thickness was 4000 A. As very often observed for 
thick YBa2Cu307_5 films, electron microscopy revealed 
the existence of well defined large inclusions, emerging 
from the flat surface of the film. Such inclusions imposed 
to increase the Ga ions FIB dose to cut the sample along 
the set of two long lines leading to the bridge itself, but 
they were not a problem for the bridge itself, which was 
patterned using two short cuts intersecting the long ones, 
in some selected area of the film free from any inclusion. 
The short cuts defining the bridge were made from one 
hundred sweeps of the FIB, etching the film down to the 
substrate (Fig. [T|). While the geometrical size for the 
bridge was about 0.2 x I — 0.8 /^m^, the actual conduct- 
ing width was likely much smaller, due to the material 
degradation along the FIB lines. The corresponding up- 
per limit for the resistivity was p{300K) < lil cm. 

Starting from the optimally doped film, the resulting 
bridge was strongly insulating. Cycling the sample tem- 
perature from the ambient down to 4.2 K resulted in 
a systematic increase of the resistance, and the results 
shown below are for one of such cycles. As can be seen 
in Fig. [21 one or the other of the 2D and 3D VRH law 
(respectively R cx exp(To/T)i/3 and R oc exp(To/T)i/4) 
may describe the data above T ~ 45 K, but the resistance 
is systematically smaller than expected for lower temper- 
ature (in what follows, we will adopt a two-dimensional 
description, owing to the increased anisotropy for YBCO 
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FIG. 3: Non linear conductance, in a representation where 
a 2D VRH gives a straight line at low electric field 
and zero temperature. From bottom to top : T = 
10, 14, 18, 26, 29, 34, 40, 46, 50, 58 K. Inset : the same data, in a 
representation similar to the one used in ref. The dashed 
lines show the determination of a characteristic electric field 
at T = 29 K. 



at low doping[18|, however stressing that both behaviors 
are indiscernible from the present data). The character- 
istic temperature in this case is Tq ~ 1.5 10"* K from the 
data in Fig. Actually, the measured resistivity was 
strongly non-linear at low temperature (Fig. [3]) and this 
might explain the bending of the data in Fig. [21 It is in- 
deed well known that the effect of an electric field in the 
VRH regime is to increase the effective temperature of 
the carriers (which determines the site occupation, using 
a Fermi distribution at the corresponding temperature), 
crossing from T^ff — T at low field to ksTef eEa 
at high field, where a is the localization lengthfl^. Our 
data is similar to the one in ref. [20| . where a complete 
numerical computation of the effect of the electric field 
was given. An estimate for the localization length may 
be obtained, using for the characteristic electric field at 
T = 29 K, E{T) = U/l ~ 1.2 10^ Vm"! (Fig[3D, and 
Ec ~ ksT/ea ([lo, [U), yielding a ~ 200 A. We remark 
that the differential resistance may appear to present a 
threshold voltage (Fig. [31 inset), as proposed in ref. p7j . 
In our case, this is a consequence of the non linearity for 
electrical field E ^ Eq. In particular, we rule out the 
Coulomb blockade as the origin of the non-linearity, as 
a tentative threshold voltage, U ~ 0.2 V (Fig[3]), would 
correspond to Coulomb energies ~ 10'^ K, whereas the 
I-V characteristic strongly evolve at much lower temper- 
ature. 

Non linear transport phenomena must be systemati- 
cally confronted to a possible contribution of Joule heat- 
ing, and one should check whether the actual tempera- 
ture of the thin film sample could be here significantly 
different from the one of the substrate. First, we shall see 
below that the localization length derived from non lin- 
earity agrees with the one derived from the dynamics of a 
single fluctuator, whereas Joule heating could potentially 
provide just any value for a. Also, fluctuator switching 
events are found thermally activated down to the lower 
temperature explored, whereas a saturation would be ex- 
pected for significant heating of the sample. Finally, we 
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may estimate the overheating, in a way similar to what 
was done in ref. p3 | . The temperature of the bridge may 
be evaluated from Tbr ^ T + 2E'^ S[p{Tbr)K]^^ , where Tbr 
is the bridge temperature, S its cross section, p the mate- 
rial resistivity and k an effective heat conductivity for the 
substrate. Using the resistance of the bridge, the heating 
is actually independent of the effective cross section. The 
high resistance helps limiting the temperature increase in 
our case. Solving this equation, using the temperature 
independent conductivity k = 150 mW/K cm and the 
resistance data in Fig. [21 we find that the heating does 
not exceed 1 K at 30 K for the higher electric field of our 
data. 

The recording of the voltage during a temperature 
sweep in the non linear regime revealed clear jumps, be- 
tween two discrete values, in a telegraph-like fashion. 
Voltage time recording, at a fixed temperature and driv- 
ing current, confirmed the existence of a voltage tele- 
graphic noise. The size of the voltage jumps were essen- 
tially temperature independent. Also, the voltage jumps 
magnitude was roughly independent from the driving 
current (Fig. [4]). However, the switching rate between 
descrete states was strongly dependent upon the volt- 
age (alt. the current) applied to the constriction, larger 
voltage yielding larger switching rate (Fig.[3]). Within ex- 
perimental uncertainties (which essentially arise from the 
difficulty to extract the telegraphic noise from the mea- 
surement noise), the switching rate is found thermally ac- 
tivated at all temperature investigated (Fig. O . We note 
that the base temperature, rather than Te//, should enter 
the Arrhenius plot for the determination of the activation 
energy, as the trapped carriers do not experience the elec- 
tric field acceleration. On the other hand, the observa- 
tion of a strong influence of the voltage on the switching 
frequency implies that the electric field work, over a dis- 
tance of the order a, lowers some activation energy for 
trapping, which is then expected of the form Eq — eaE. 
Using the data in Fig. [H we obtain from the slope of the 
linear fit of the inset, dLn{f)/dU = ea/{l ksT), a = 170 
A and Eo/kB = 220 K, the localization length being in 
full agreement with our former independent determina- 
tion. It should be stressed that we did not observe any 
switching event in the linear regime evidenced in Fig.[3l if 
such switching events exist in the linear regime, the rate 
must be so low for the currents involved for this regime, 
that it is out of our experimental time scale. As a con- 
sequence, we believe that the telegraphic noise observed 
here is characteristic of the non linear regime, and can- 
not easily be interpreted using standard linear models for 
random telegraph noise |23i although the origin is likely 
the same : the capture of carriers into localized defects 
which either directly remove carriers from the conduc- 
tion band, or alter the mobility of the carriers. Many 
such events, in larger samples, may be at the origin for 
1 // noise in the vicinity of the superconductor-insulating 
transition Q. 

For some recordings, the signal exhibited three states 
histograms, as shown in Fig. [S) The simplest interpre- 
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FIG. 4: Mean voltage jump for different driving voltage. In- 
set: switching rate, showing thermally activated behavior. 
The line is a linear fit, yielding a = 170 A(r = 25 K.). Note 
that the switching rate becomes unmeasurably small upon 
entering the linear regime displayed in Fig. |3l 
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FIG. 5: Arrhenius plot of the switching frequency {1 = 1 ^A). 
The line is the best fit to thermally activated behavior, with 
activation energy 13 K. 



tation is that such a signal results from the presence of 
two distinct traps, which independently capture an re- 
lease charges. However, we systematically observed that 
the signal level separation was nearly the same between 
adjacent levels. This could have been only fortitious, 
but we have observed that the levels occupancy evolved 
with driving current in a correlated manner: for inde- 
pendent traps, such an increase should simply increase 
the switching rates of both traps, whereas we observe 
that the weight of the distribution is progressively trans- 
fered from the lower level to the upper one (Fig. [S]). This 
correlated behavior suggests that, rather than two inde- 
pendent traps, one actually has two traps correlated in 
some way, or a single trap with three states. 

Such a case was first reported for the linear conduc- 
tance of a MOSFET channel[24]. As underlined by the 
authors, the correlated two traps case requires that these 
should reside at a close distance, so that Coulombic in- 
teraction between them can shift the energy level of the 
trapped charge by several fc^T : this implies cluster- 
ing of the defects at point-like nucleation sites, at a dis- 
tance that does not, however, allow for hybridization of 
the defects. The second possibility is the existence of 
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FIG. 6: Left : three states voltage distribution for driving 
current / — 0.7, 0.8, 0.9 fj,A (the origin for U was arbitrarily 
shifted) . The weight of the distribution is progressively trans- 
fered to the upper level with increasing current. Right : part 
of the time recording used to establish the distribution. 

traps which can accommodate multiple charges, as re- 
ported in MOSFETsfii, [H- The nature of such three 
states traps is however unclear. In the case of Anderson- 
localized states, the possibility for doubly-occupied local- 
ized states exists as, when the on-site correlation energy, 
U, is smaller than the bandwidth of the localized states 
distribution, some of the localized sites may be doubly 
occupied, with two statcs[26]. The upper state energy is 
close to Ep, as for singly occupied states, while the lower 
state energy is about U below the Fermi level. How- 
ever, for the higher electrical field reached in the present 
study, E/Eo ~ 0.1, we have Tefj/To ~ 0.1 ~ 10^ K, 



which is too low to allow for the activation of the state 
lying typicallyjl^ ~ 2 eV ~ 2 10^ K below Ep. Then, 
transport in our moderate out-of-equilibrium conditions 
should involve only the higher states, as is the case for 
the linear transport regime ^2^. Finally, one should also 
consider the possibility for a contribution of collective 
charge motion. Indeed, it has been proposed that charges 
may gather into conductive channels - or '"stripes'" - as 
a Mott insulator is doped. As underlined in Ref. [23 |. 
collective motion of charges generally implies a sharp 
threshold electric field and a narrow band noise. The 
latter property should however be characteristic of con- 
ventional charge order extending along large coherence 
length, rather than of quasi-lD stripes, as expected in 
cuprates, so that what one really would expect, in the 
present case, are a low threshold electric field and large 
conductance fluctuations, associated to the existence of 
a large number of correlated charges. These two charac- 
teristics were both observed in the present study, but we 
have favored a more conventional interpretation based on 
transport in the VRH regime. We find that it is able to 
account for several of our observations but, clearly, the 
exact mechanism for some correlation that we have ob- 
served for the trapped charges needs to be investigated in 
a more systematic way. While our tuning parameters for 
the dynamics of an individual fluctuator was restricted 
here to the voltage that tunes the accessible states in the 
non-linear VRH regime, and the temperature that allows 
for the relaxation of the charges from traps, it would be 
desirable to be able to investigate the behavior of indi- 
vidual fluctuators in the linear regime also. Reversibly 
tuning the doping of the material over some short range, 
using the field effect on very thin films, just as one tunes 
the channel conductivity of a MOSFET from the gate 
voltage, appears as a convenient way to do so. 
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